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Abstract 

Inhomogeneity is introduced through random local interactions 
(Ui) in an attractive Hubbard model on a square lattice and studied 
using mean-field Bogoliubov-de Gennes formalism. Superconductivity 
is found to get suppressed by the random Ui contrary to the results of 
a bimodal distribution of Ui. The proximity effect of superconductiv- 
ity is found to be strong, all sites develop non-zero pairing amplitude. 
The gap in the density of states is always non-zero and does not vanish 
even for strong disorder. When two such superconductors are coupled 
via a channel, the effect of one on the other is negligible. The length 
and width of the connector, do not seem to have any noticeable ef- 
fect on the superconductivity in either systems. The superconducting 
blocks behave as independent entity and the introduction of the chan- 
nel have no effect on them. 

PACS: 74.81.-g, 71.10.Fd, 74.20.-z 

1 Introduction 

One of the most important aspects in the study of correlated system is the 
study of the effect of spatial inhomogeneity p. There have recently been 
several examples of such systems where these inhomogeneities can occur in- 
trinsically via quenched disorder in the system, or can occur spontaneously. 
For example, impurities can be driven in a superconductor by irradiation 
or chemical substitution. On the other hand, holes in the cuprate super- 
conductors or magnetic or charge ordered domains in manganites sponta- 
neously arrange in geometric patterns like stripes or a checkerboard at cer- 
tain fillings [3, SI El 13 . Using scanning tunneling microscopy evidence for 
electronic inhomogeneity has been reported in the high-Tc superconductor 
Bi2Sr2CaCu20s+x by McElroy, et al. [2j. Most of the high-Tc materials are 
ceramic in nature and inhomogeneities are present in even the best prepared 
samples. This inhomogeneity is manifested as spatial variations in both the 
local density of states and the superconducting energy gap[7j. 

There has been a number of theoretical attempts to understand the effect 
of quenched disorder in superconductors P, ^ . In the context of repulsive 
models, both in the weak-coupling models and their strong-coupling coun- 
terparts, considerable numerical work has been done with inhomogeneities 
p!Ol \Tl\ . In such cases stripes and checkerboard patterns have been reported 



[I2l[l3|, the presence of d-wave superconductivity, however, is less convinc- 
ing. Enhancement of Tc due to inhomogeneities in the weak-couphng regime 
is demonstrated by Martin et al. [14j. Aryanpour et al. [8] studied an s-wave 
superconductor with quenched disorder starting from a negative-U Hubbard 
model using a mean-field theory. The disorder enters in their model through 
a random choice of two values of the attractive interaction (bimodal distribu- 
tion) at different sites. Quite interestingly, it was shown that below a certain 
value of the average attraction, the zero temperature superconducting gap 
is larger than that of the homogeneous superconductor with same (uniform) 
attraction. Ghosal et al., [9j use an s-wave superconductor and look for 
the eflFect of disorder using a mean-field treatment. Both these calculation 
observe strong effects of disorder on the superconducting order parameters. 
In certain specific cases, lower dimensional orders, like stripes or checker- 
board order [5l [15] has been found [8j in the numerical calculations. One 
possible application of these inhomogeneities is to manipulate the degree of 
superconductivity and the transition temperature through varying degree of 
inhomogeneities. Another possibility is the study of the coupling of two ad- 
jacent systems, where superconductivity in either could be controlled by the 
degree of disorder in them. In many of the tunneling devices such conditions 
are obtained. 

We consider a similar situation with the local attractive interaction taken 
completely random, which, we believe, is more realistic than using a bimodal 
distribution. We take a system of two such disordered, superconducting 
blocks connected by a channel. This model conforms to various actual sit- 
uations in tunnel junctions and the study of them is quite interesting. We 
look into the problem of disordered superconductivity and compare the re- 
sults with homogeneous systems. Using mean-field exact diagonalization 
techniques, we work out the case of random disorder (as opposed to forcing 
specific geometrical patterns like stripes [8j from the outset) in a two dimen- 
sional square lattice and also for a system of two such blocks with a channel 
connecting them. We compute various microscopic parameters such as the 
local pairing amplitude A^, the average electron density n^, the local chemical 
potential /x^ and the density of states for different fillings for various average 
interactions. We also look for any emergent one dimensional pattern in the 
real-space. The effect of coupling of two such systems is studied in detail. 

2 Model and calculation 

Our starting point is the attractive Hubbard model given by 

H = -t i^cjCja + c]^Ci^) - /x^cJ^Q^ - Y\Ui\ni^nii (1) 

<ij>a i i 

where t is the hopping potential, /x is the chemical potential and Ui is the 
local attractive interaction between the fermion of the opposite spins residing 



on the same lattice site i. c]^ and Cia- are creation and destruction operators 
for an electron with spin a on site z, hi^ and are number operators at site 
i with spin up and down. Superconductivity and other charge instabilities in 
the homogeneous case have been studied in great detail [161 EZl the past. 
At half filling the superconducting state is degenerate with a charge density 
wave state and they are connected by a pseudospin rotation. 

In order to study the effect of disorder, an inhomogeneous calculation is 
necessary. We use the Bogoliubov-de Gennes (BdG) mean-field approxima- 
tion and replace local electron correlation by the local superconducting pair- 
ing amplitude A^, at site i (only s-wave pairing considered here). = {ci^Cii) 
and = cl^Cia. Assuming {ni^) = {nii) = {ni)/2we get, 



^eff = Yl i^aCja + c]^Ci^) - fiiclc,^ ^ |[/^| [A.cJ^cJ^ + A*C,|Q^] (2) 
<ij>,cr i(7 i 

where fii = /a + |[/^|(n^)/2 is a site-dependent Hartree shift with (n^) = 
Zlcr i^ia)- All energies are scaled to t = 1. We use an initial value of the 
chemical potential fi and use different random configurations of Ui to realize 
the inhomogeneous superconductor. We then compare its tendency of super- 
conductivity with the homogeneous system, that is, Ui = U for all i. For the 
inhomogeneous case we have taken Ui to be uniformly random between two 
specified values. In all these cases we use the self-consistent diagonalization 
of the Hamiltonian - at every iteration compute the averages at every site, 
recalculate the Hamiltonian and set out for diagonalization again. When 
the iterations converge we obtain a self-consistent result as usual and then 
compute the values of the physical the quantities. 



3 Results and discussion 

We take a square lattice of finite size and set up the BdG Hamiltonian. On 
reaching self-consistency, the average electron density, i.e., the average filling, 
is calculated. This is repeated for many initial fi and the results are inverted 
to find the necessary physical quantities for a fixed (n) as is customary in 
grand canonical ensemble (the local inhomogeneity and self-consistent inter- 
actions forbid fixing the average electron density). We note that ideally one 
would expect only topological order in two dimensions. However, this does 
not apply to mean-field order. 

The variation of A^ with Ui is shown in Fig.l. As expected, the sites 
with large negative attraction have larger superconducting amplitude. It is 
interesting to note that although there are sites with strictly zero or very 
small \Ui\^ there are no sites with zero A^. This is not unexpected, as a 
large local fiuctuation of order parameter is quite unfavourable with respect 
to quantum fiuctuations and are energetically costly in terms of the loss of 
interaction energy. Such finding have also been reported earlier [8j, but we 



do not find any value of {U) for which the disordered system has a larger 
average gap than the uniform case. In the case of a bimodal distribution of 
Ui^ though, the proximity effect seems to be much stronger and the disordered 
system has a larger gap than the uniform one below a certain ([/). We also 
observed that unless a specific geometry of the disorder is quenched into the 
lattice, there is no possibility of a stripe like state to organize self-consistently 
as observed by DagottofSj. We believe in a situation where the disorder is 
annealed and randomly distributed, such highly anisotropic states are quite 
unlikely in the present model. 




Figure 1: (Colour online) (a) Variation of the local pairing amplitude with 
the interaction [/^, distributed uniformly in the range 0-4 (b) Comparison 
between the variation of < A > with < [/ > for uniform and disordered 
cases at < n >= 1, < n >= 0.8 and < n >= 0.6. (c)Variation of < 
^uniform > / < ^disorder >with < U >. (d) Variation of average local pairing 
amplitude< A > with < n >)(e) Superconducting regions for U=0-2. The 
red islands represent superconducting regions with large amplitude (see text). 



Fig. la shows that increases with Ui and saturates at large [/^, there is 
a broad distribution of for lower values of Ui. Fluctuation in the pairing 
amplitude A^ are indeed larger in the weak-coupling region. The variation of 
the average Ui and the average < A >= j^Y^i gives a good indication of 
their relationship and can be compared with the uniform case. The compar- 
ison in Fig. lb shows an enhancement of the average superconducting pairing 
amplitude (A) with (U) as expected saturating at a higher value of the at- 
traction. The value of average pairing amplitude is higher for the uniform 
case than the disorder case for all values of< n >= 1.0, 0.8 and 0.6 contrary 
to that of bimodal distribution. Fig. Ic shows the ratio as a function of av- 
erage interaction, and it rises sharply at low U and always stays above one. 
This implies that the random disorder inhibits long range superconducting 
order. Fig. Id shows the usual bell-shaped curve of < A > versus < n > with 
maximum at < n >= 1. 

In order to glean the real-space picture of the superconducting regions, 
we mark the sites where the value of A^ is greater than the corresponding 
value in the homogeneous case (Fig.le). Since all sites have non-zero A^, this 
prescription allows us to locate preferred patterns, if any, in the real space. 
We do not see any such pattern while the regions of strong A^ are quite ran- 
domly distributed forming patches. The formation of islands in turn reduces 
the overall observed Tc (which is determined by the phase coupling across 
such regions) in a real system and the corresponding superfluid stiffness. 

Plotting the frequency distribution of A^ gives a good account of the ef- 
fects of disorder. The peak in Fig. 2a is certainly above the mean value. 
With higher disorder the peak shifts towards right. Further on we calculate 
the pair correlation function (Fig. 2b), defined as c(r^ — rj) =< A^ Aj >. 
A disorder averaging is done as usual to restore the translational invariance. 
The correlation function drops sharply at short distances but saturates at 
large distances indicating a true (mean-field) long range order. The correla- 
tion function for four ranges of U are shown in Fig. 2b, they all merge to a 
single curve (except for small fiuct nations due to finite size) when normalized. 




Figure 2: (a) Frequency Distribution P(A) of A for U=0-2 (b)Plot of corre- 
lation for four U ranges: 0-2 0-4 0-6 0-8 
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Figure 3: Density of States plot for U=0-2 and 0-4, the two sets of curves 
are offset vertically to make them separate. 

The presence of a superconducting gap is another indication of the true 
long range order [9j. It is, therefore, very important to look at the density 
of states (DOS) of the system. As mentioned above there are no sites with 
zero pairing amplitude. This is reflected in the DOS where the gap at zero 
energy is clearly seen in Fig. 3, while the very sharp peaks (divergences in 
a homogeneous superconductor) are now broadened into two symmetrical 
regions with a broad distribution of states at higher and lower energies. The 
presence of the energy gap even in the highly disordered systems clearly 
indicates that s-wave superconductivity is possible even in the presence of 
large disorder. 

To understand the effect of coupling of two disordered superconductors, 
we take up two such superconductors and join them via a narrow channel 
(Fig. 4f) and observe the effects due to the proximity of each other. We 
find that contrary to the common expectation the width and length of the 
channel do not have any significant effect on the superconductivity in either 
systems (4d,e).More interestingly, results from the comparative variation of 
(A), ([/), (n) of the two blocks (Fig. 4a-e) do not give any clear infiuence of 
the channel. For example, we observed the variation of (A) of one supercon- 
ductor with the {U) in both uniform and disordered case and compared with 
similar results of the single block system we had obtained earlier. Results are 
symmetric with respect to interchange of systems 1 and 2. It is clearly seen 
from Fig. 4a that the average superconducting pairing amplitude in blocks 
of the coupled system follow identically the pattern of a single block system 
in both disorder and uniform regime. The channel was kept homogeneous. 
The results do not change if the channel is maintained homogeneous but non 
superconducting {Ui = for all i in the channel). This unexpected behaviour, 
we believe, is due to the formation of islands of low pairing amplitude dis- 
cussed earlier (Fig le). These regions localize the electron in the block since 
it would be energetically unfavourable for the electrons to percolate through 
the channel. This inhibits correlations between the blocks, effectively making 



them behave as independent systems. 

The average superconducting pairing amphtude in the second system in- 
creases at the expense of (Ai) (Fig. 4b since the enhanced attraction in 
system 2 provides extra stabihzation energy in that region. Note that the 
variation of < A2 > with < C/2 > is also unaffected (Fig. 4b) by the presence 
of system 1 and nearly similar to that for a single system. The variations 
of overall average pairing amplitude (Ai + A2) with total average interac- 
tion at a fixed overall density is shown in Fig. 4c. The combined system 
behaves like an isolated single system and expectedly shows the typical rise 
and saturation behaviour seen in a single system. 

Fig. 4d demonstrates clearly how the average order parameter in the two 
systems are correlated with the corresponding interactions. The larger the 
attraction, more is the average (A). As (U) increases in any system, the 
density of electrons tends to increase for the extra stabilization available in 
that region. Fig. 4e shows this tendency, with (n^) in each superconductor 
increasing rapidly with the average attractive interaction there. 

4 Conclusion 

In conclusion, we have studied inhomogeneous s-wave superconductors us- 
ing the Bogoliubov de-Gennes mean-field theory. Superconductivity is sup- 
pressed over the homogeneous weak-coupling value due to disorder, though 
the proximity effect is strong with non-zero pairing amplitude at all sites 
(even with Ui = 0). The frequency distribution ofp{A) shows a peak towards 
higher values of A. The gap in the density of states persists even for the high 
disordered case lending support to the existence of strongly disordered in- 
homogeneous superconductivity. When two such blocks of inhomogeneous 
superconductors are connected via a channel, there is no appreciable effect 
of one on the other. The individual blocks are not affected significantly by the 
other block connected by the channel . We argue that this is due to the fact 
that the fiuctuating local order parameters remain pinned to the individual 
values thereby preventing significant correlation between the blocks. 
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Figure 4: (a)Variation of < A > with < U > for individual blocks of 
the coupled system for both disorder and uniform case and its comparison 
with single block system, (b) (A)2 as a function of {11)2 for {71)2 = 1 and 
comparison with superconductivity in a single system, (c) (A) as a function 
of (U) with (n) = 1 for the whole system, (d) Variation of F = (A)2 — (A)i 
with q = {11)2 — {U)i^ {n) = 1. Here the legend / refers to length of channel 
and y refers to its width, (e) Variation oi u = {n)2 — {n)i with q while 
(n\ — ^. (f) Two snnermnductinp' blocks connected bv a, channel 



